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Abstract 

We present a new model for the propagation of polarized light in a random birefringent medium. 
This model is based on a decomposition of the higher order statistics of the reduced Stokes parameters 
along the irreducible representations of the rotation group. We show how this model allows a detailed 
description of the propagation, giving analytical expressions for the probability densities of the 
Mueller matrix and the Stokes vector throughout the propagation. It also allows an exact description 
of the evolution of averaged quantities, such as the degree of polarization. We will also discuss 
how this model allows a generalization of the concepts of reduced Stokes parameters and degree of 
polarization to higher order statistics. We give some notes on how it can be extended to more general 
random media. 

1 Introduction 

The interaction of polarized light with a random medium is of interest to many domains. Examples 
include imagery, telecommunications, medicine and instrumentation. In this paper, we present a new 
model for the propagation of a polarized lightwave in a random birefringent medium. This problem 
is relevant to telecommunications. Indeed, it occurs in optical fibers subject to Polarization Mode 
Dispersion, or PMD. See for instance [2,8,17]. The review [11] gives a comprehensive introduction to 
PMD. The approach offered in [2, 17] is based on stochastic differential equations satisfied by the Jones 
(or Mueller) matrix of a certain length of the medium. 

Our main aim in this paper, is to introduce a new model for the propagation of polarized light in 
a birefringent random medium. All the main features of this model are generalizable to random media 
in general. However, we use this simpler situation as a starting point. Unlike the ones given in [2, 17], 
the model we present is global and not local. This means that instead of using a stochastic differential 
equation that describes propagation over short "infinitesimal" distances, we concentrate on the linear 
operators describing propagation over finite distances. The two approahces are complementary. However, 
the results that we will describe are more far-reaching than the ones obtained using stochastic diffenretial 
equations. 

One of the main features of this model is that it treats the random medium as a whole. In other words, 
the random medium is considered as a system and a certain type of statistical relation between its input 
and output is studied. The microscopic characteristics of the medium are thus related to the parameters 
of the probability laws involved in the input-output relation characterizing the medium (system). They 
can be estimated given samples of the input and output of the medium. In fact, the whole probability 
density of the Mueller matrix of the medium can, in principle, be estimated given samples of the input 
and output. 

The model we will present is based on a decomposition of the higher order statistics of the reduced 
Stokes parameters along the irreducible representations of the rotation group. The rotation group is 
important because a birefringent medium acts on the Poincare sphere by rotations [3]. The relation 
between the probability densities, on the Poincare sphere, of the input and output reduced Stokes 
vectors will be shown to be given by a spherical convolution. The decomposition along irreducible 
representations of the rotation group is used to obtain from this relation the laws of transformation of 
the higher order statistics of the reduced Stokes vectors. On the one hand, these laws of transformation 
constitute generalizations of the Stokes-Mueller formalism to higher order statistics. On the other hand, 
they will be shown to be a powerful tool for describing the propagation of polarized light in a random 
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medium. This whole approach is actuaUy an example of harmonic analysis on the rotation group [13] 
which is an instance of non commutative harmonic analysis [10, 13,25]. 

Classical models for the state of polarization of a lightwave, such as the Jones and Stokes models, 
only involve statistics of order 1 and 2 [3] . Whenever the fluctuations of the field of the lightwave are non 
Gaussian, it is necessary to consider higher order statistics. Non Gaussian fields occur in many situations 
in optics [5]. Different ways of including higher order statistics have been proposed. For instance, 
Refregier in [6, 24] gives a measure of the degree of polarization based on KuUback relative entropy 
between the field of the lightwave and a completely isotropic field with the same intensity distribution. 
In [16], Luis proposes a measure of the degree of polarization based on the mean quadratic distance 
between the probability density function on the Poincare sphere of a lightwave and the probability 
density function of a uniform distribution on the Poincare sphere. In [5] Ellis and Dogariu propose using 
the correlations of the Stokes parameters to discriminate states of polarization with the same second 
order statistics. In [4], they use the symmetries of the probability distribution of the Stokes parameters 
on the Poincare sphere to make this same distinction. All these models emphasize two aspects. That 
it is important to include higher order statistics and that this should be done in accordance with the 
Symmetries of the probability distribution on the Poincare sphere. Our approach is quite similar to them. 
It has the additional advantage of formalizing the importance of the spherical symmetry underlying the 
problem, by using the irreducible representation of the rotation group. 

The plan of the paper is the following: In section [2] we outline the physical situation we wish to 
consider. In section [3l we explain the mathematical tools we wish to use, namely the irreducible rep- 
resentations of the rotation group. In section [4] we give the main equations of our model. In section O 
we discuss how it generalizes the Stokes- Mueller formalism to higher order statistics. In section [5] we 
show how it can be used to describe the propagation of a lightwave in a random birefringent medium. 
In section [71 we discuss how our model can be extended to the case of a general random medium and 
conclude. 

2 Physical problem 

Let us quickly outline the kind of physical setup we are interested in. We consider the effect of a random 
birefringent medium on the state of polarization of an incident lightwave. We will be interested in two 
kinds of relations. The first kind of relations is between the input and output states of polarization when 
the medium is taken as a whole. The second kind is the relations governing the changes in the state of 
polarization of the lightwave during its propagation in the random medium. The second kind of relations 
can be considered as a more detailed model for the first. 

In general, the state of polarization of a lightwave will be modelled as a random variable on the 
Poincare sphere [3] . If we are considering the input and output states of polarization to a birefringent 
random medium taken as a whole, we will note Sin and Sout the random variables on the Poincare 
sphere corresponding to the input and output states of polarization. The effect of the random medium 
on Sin is given by its Mueller matrix. Here, it is understood that this Mueller matrix is a matrix-valued 
random variable [26] . In general, the Mueller matrix acts on the Stokes vector of the incident lightwave 
and not on the reduced Stokes vector Sin- However, in the special case of a birefringent medium the 
Mueller matrix reduces to a rotation matrix acting on the reduced Stokes vector Sin ■ According to this 
discussion, the effect of a random birefringent meduim is given by a random variable R with values in 
the rotation group 50(3). This effect is given by the following equation: 

Sout — RSin (1) 

Where Sin and Sout are random variables on the Poincare sphere S^ and i? is a S'0(3)-valued random 
variable. 

When considering the propagation of a lightwave in a birefringent random medium we will give its 
state of polarization by a stochastic process S{z) with values on 5^. Here z > is the distance along 
which the wave has propagated through the random medium. For every z, S{z) is a random variable 
on S'^ giving the state of polarization of the wave after a distance z of its trajectory in the medium. In 
relation to equation H]), we can note S{0) = Sin and S{Z) = Sout for some given distance Z. We make 
the hypothesis that a length z of any trajectory of the wave can be described using a rotation matrix 
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(as in II])) noted R{z). We have, as a first model for the propagation: 

S{z) = R{z)S{0) (2) 

By associating a matrix R{z) to every distance z, over which the wave propagates in the medium, we are 
imphcitly ignoring any scattering or beam divergence effects in the medium. This amounts to supposing 
that hght remains collimated in the medium [18]. 

In the following section we give the mathematical tools we will use to analyze equations ^ and ^ 
to obtain the main equations of our model as given in section |H 



3 Irreducible representations of the rotation group 

The main idea of this paper is to decompose the higher order statistics of the reduced Stokes parameters 
along the irreducible representations of 5*0(3). This decomposition is central to the spherical convolution 
theorem which we present shortly and which allows us to obtain the laws of transformation of the higher 
order statistics of the reduced Stokes parameters. These laws of transformation constitute the main 
equations of our model and are given in section ID In fact, we will be interested in the realization of the 
irreducible representations of 5*0(3) in the space of square integrable functions on S^. Under this form, 
the study of the irreducible representations of 50(3) is known as spherical harmonic analysis [10, 13,25j3- 
Let /(s) = f{(j), 6) - where (j) and 9 are the azimuth and polar angle of s = (si, S2, S3) - be be a square 
integrable function on 5*^. / can be decomposed along the orthogonal basis of spherical harmonics Yj^{s), 
where I € N and —l<m<l: 

m—l 

/W = E E i2l + ^)fLYi{s) (3) 

Where /4, is the projection: 



fL = / fis)Yiis)ds (4) 

ds is the Haar measure ds — sin <j)d(j)dd and the bar represents complex conjugation. Remeber that 
spherical harmonics are given by the formula: 



Where the are associated Legendre functions. And also by the cartesian formula: 



^^mV^E^I±^ (^)"/-(.3+^yr^coso'cos(™i)di 

Under the effect of a rotation r E 50(3), the (2/ + 1) coefficients fl^ (for every I e N) transform according 
to an irreducible unitary representation of dimension (21 + 1) of 5*0(3). In other words, for the rotated 
function /r(s) — firs) we have the development: 

m—l 



Ms)^Yl E (2; + l)(/.)„i^,Us) (7) 

;eN m=-i 

» I 

Where the coefficients {fr)m of the development are given (for every Z e N) by the linear transformation: 

I n—l 

ifr)m = E ^DOO/^ (8) 



^See these same references for all the results in this section. 



3 



The functions on SO{3), D\^^^{r), which give the elements of the matrix of the linear transformation ([5]) 
are the matrix elements of the irreducible representation of dimension {21 + 1) of SO{?>). They can be 
given explicitely, as functions of the Euler angles 6*, ?/;) of the rotation r: 

DLn (r) = {<t>,e,rP)^ e-™^Pi„ (cos e)e~^-* (9) 

Where the polynomials are given by Jacobi polynomials. 

The {21 + 1) X {21 + 1) matrices Z)'(r) = —l<m,n<l, are unitary and verify the 

homomorphism property: 

D'{nr2) = D'{n)D'{r2) (10) 

Where ri,r2 G 50(3) and I is any natural number. 

The functions Dl^^{r), on SO{3), play a similar role to that played by spherical harmonics on S^. 
Indeed, let g{r) be a square integrable function on SO{3), then g can be decomposed along the orthogonal 
basis of matrix elements D\^^{r): 

m.n—l 

5W=E E i'^l+^WmnDLJ^) (11) 

iSN 771,11 = — I 

Where gj„„ is the projection: 

gLn = / f{r)Dl^{r)dr (12) 

JS0{3) 

and dr is the Haar measure on 5*0(3), given in terms of Euler angles by dr = snv6d(j)d0dij). 

It is interesting, at this point, to rewrite the developments ([SJ and (jlip in a matrix form. This will 
clarify the rest of the article and avoid cumbersome notations. To do this, adopt the following notations: 
Note for every / G N, the {21 + 1) dimensional column vector whose components are the coefficients 
fmy~l^''^^h from equation Q. In the same way, note Y^s) the column vector of spherical harmonics 
YJ^{s). Note the {21 + 1) x {21 + 1) matrix with elements gj„n; equation ([TT|) . The notation D^{r) 
has allready been explained. Using these notations, equation ^ can be rewritten as follows: 

/(.) = ^(2? + i)(/')*r'(s) (13) 

Where t stands for transposition. Equation pip can be rewrittin as: 

g{r)^J2(^l + l)Tr[g'{D'{r)y) (14) 

Where Tr stands for the trace and f for the Hermitian conjugate. 

As explained above, the spherical convolution theorem will be used to obtain the basic equations of 
our model. This theorem uses the decompositions ()13p and ()14p to transform a spherical convolution 
into a family of matrix products. It is an analogue of the classical convolution theorem, which is used to 
transform a convolution of functions of a real variable into a product of their Fourier transforms. 

The convolution of two functions 52 and 171 on 50(3) is a function g on 50(3) defined as: 

gir) = (52 * gi){r) = f g2{t)gi{t-\)dt (15) 

J 30(3) 

This definition is analogous to that of the convolution of two functions of a real variable. Formally, it 
can be obtained from it by replacing the usual r — t by t^^r. The convolution of a function g on 50(3) 
with a function / on 5^ is a function h on . It has a similar definition: 

h{s)^{g^f){s)^ I g{t)f{t-\s)dt (16) 

J 30(3,) 

The convolution theorem states that equation ()15p is equivalent to: 

g' = gig[ for I e N (17) 
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Where the product on the right hand side is a product of {21 + 1) x {21 + 1) matrices. It also states that 
equation (jl6p is equivalent to: 

= gip for / G N (18) 

Where the product on the right hand side is a product of the {21 + 1) x {21 + 1) matrix g' with the (2^ + 1) 
dimensional vector 

4 Main equations of the model 

In this section, we give the main equations of our model. In the following section [5l they will be 
interpreted as giving the laws of transformation of the higher order statistics of the reduced Stokes 
parameters. In section [51 they will be used to give a detailed probabilistic description of the propagation 
of polarized light in a random birefringent medium. 

Let us return to the physical situation described by equation ([T]) of section [21 This equation describes 
a random medium taken as a whole. It gives the relation between the random variables on the Poincare 
sphere, Sin and Sout , describing (respectively) the input and output states of polarization to the medium. 
In this section, the mathematical tools given in the last section [3] will be used to analyze equation ([1]). 
First, We will establish the relation, between the probability densities of the random variables Sout, R 
and Sin appearing in this equation, as a spherical convolution. Then we will use the convolution theorem 
([T8|) to transform this relation. 

Note Pin{s) and Pout{s) the probability densities of Sin and Sout- Note PR{r) the probability density 
of R. To see that the relation between these three probability densities is given by a spherical convo- 
lution ((TH)) . apply the law of total probabilities: The probability for Sout to take a value near s G S'^, 
conditionally to the event " R is near r e 5*0(3)" is equal to pin{r~^s). Indeed, the random variables Sin 
and R are independent (the properties of the medium do not depend on the incident wave). By summing 
these condition probabilities, we obtain: 

Pout{s) ^ / PR{r)pin{r^^s)dt = {pR* pin){s) (19) 
JsO(3) 

Which is a spherical convolution. This relation is analogue to the classical relation stating that the 
probability density of the sum of two independent real random variables is the convolution of their 
probability densities [12]. 

We now apply the convolution theorem - equation (jlSp - to this last relation. Note and ^^^j the 
coefficient vectors in the development ([T3| of (respectively) pin and Pout- Note R^ the coefficient matrices 
in the development of pr. Then according to P^ : 

cLt = R'^iL (20) 

This is the first main equation of our model. First of all, it describes the state of polarization of a 
lightwave using a family of coefficients, e.g. £,outJ instead of a probability density function on S'^. 
It also describes the law of transformation of these coefficients, when the lightwave passes through a 
random birefringent medium. This description presents some similarities to the model suggested in [16], 
for the desciption of quantum polarization effects. 

Now consider that Sin is affected by the composition of two "random elements". That is, let Sout — 
R2RiSin. Where R2 and i?i are, as in ([T]), random variables in 5*0(3). Let R = R2R1 and note pR{r) 
the probability density of R and similarly note pr-^ and pR^ . Applying the same reasoning as the one 
used to obtain equation (|19|) . we have: 

PR{r) = / PR^{t)pR,{r^r)dt = {pR,_ *PR,){r) (21) 
Jso{3) 

In order to apply the convolution theorem (|17p to this equation, note i?', i?2 and R[ the coefficient 
matrices in the development ()14p of (respectively) pR, pR^ and pr^^. Then, according to p7|) : 

R^ = R^^R[ (22) 
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This relation can be generalized to the composition of n random elements, R — . . . i?2^i- Using 

a similar notation to the one in I22[ we have: 

This is the second main equation of our model. It shows how equation (|20p can be applied when the 
input wave is subjected to the consecutive effect of several random elements. This relation will be used 
in section [S] where we study the propagation of polarized light in random birefrigent media. 



5 Statistical interpretation and generalization of Stokes formal- 
ism 

In this section, we give a statistical interpretation of the results of the last section, especially equation 
(PU)) . We study the relation of these results to the Classical Stokes-Mueller formalism [3] as well as to 
recent works that study the role of higher order statistics in polarization optics [4, 5, 16, 24]. The main 
idea is to notice that the coefficient vectors and ^^^j, of equation ([20| . contain combinations of the 
moments of order I of the random vectors Sin and Sout , that transform under irreducible representations 
of 50(3). This observation is used to generalize the notions of reduced Stokes vector and of degree of 
polarization to higher order statistics. 



5.1 Generalized reduced Stokes vectors 

Let us consider a random variable S on the Poincare sphere, describing the state of polarization of a 
lightwave. Let p{s) be the probability density of S. If ^' are the coefficient vectors of the development 
(fT3| of p, then by definition -see equation jl]): 



e=/ p{s)Y\s)ds = ¥.{Yi{S)) (24) 

By considering the Cartesian expression ^ for the spherical harmonics F'(s), it is possible to see that 
the coefficient vectors contain complex combinations of the moments of order I of the vector S. Let 
us take the example oil = 1. Using formula ^ we can see that the vector Y^{S) is given bjH Y^{S) — 

(^Si^^g^^_Sr^^\ It results from this that ^1 = (^^^^^^,£(53),-^^^^^)*. In other words, 

is related by a complex change of basis to the average reduced Stokes vector E(5) = (E(5i), E(S'2), E(53)). 

The degree of polarization is classicaly defined using the average reduced Stokes vector. That is, 
using the second order statistics of the field of the lightwave. It is given by ^"^^ 



P = mS)\\ = v/E2(5i) + &{S2) + &{S^) (25) 



Note that we can also write P = \\(}\\ = y IC-iP + ICoP + 

This example shows that using the first coefficient vector of the development of the probability 
density of S, we retreive the classical average reduced Stokes vector as well as the classical notion of 
degree of polarization [3]. To generahze this result to higher order statistics, we construct, for every 
Z £ N, a real version of the [21 + 1) dimensional coefficient vector . This can be done as follows. Define 
for every I the (2^ + 1) dimensional real vector 5*', as follows: 

„+eL)=-V23?(eL) for™>0 
5" = <( 5^ = for m = (26) 

'SL - d„J = %/23(e™) for m < 




^Y\s) has been defined as a column vector, whence the transpose. 

^For the definition I I25I I to correspond precisely to the usual definition of the degree of polarization, we must add the 
hypothesis that the total intensity of the lightwave is independent from the other three components of the Stokes vector. 
This is not a very restrictive hypothesis when considering birefringent media. 
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We will call the vector S the reduced Stokes vector of order I. It contains the moments of order I of the 
vector S. For example: 



(27) 



and, 



E{S). 



The definition of the vectors S allows the generalization of the notion of degree of polarization to 
higher order statistics. Indeed, by analogy with formula (|25p . we can define: 



We will call P' the degree of polarization of order I. It is possible to prove that P' G [0, 1]. Indeed, P' 
is evidently positive. Note also that: 



50 that P' < 1. For I = 1, the fact that P^ = P G [0, 1] is well established [3]. It means that the only 
physical states of polarization are the ones with E(S') inside the Poincare sphere. 

5.2 Examples and relation to other works 

In this section, we give a few examples of how the notions of reduced Stokes vector of order / and degree 
of polarization of order I, introduced in the last subsection 15.11 can be used to distinguish states of 
polarization which are indistinguishable in the framework of classical models for polarization [3]. We 
also explain how our model is related to other recent works on higher order statistics in polarization 
optics [4,5,16,24]. 

Remember that -see the introduction- classical models for polarization only use the first and second 
order statistics of the field of the lightwave [3] . This corresponds to using the average Stokes vector or 
reduced Stokes vector. This approach is sufficient for Gaussian fields but fails for non Gaussian fields [4,5]. 
Indeed, considering only the average reduced Stokes vector E{S) = would lead to identifying states 
of polarization which have the same average E(>5') but might have different higher order moments of this 
vector. That is, in the formalism introduced in the last ssubsection 15.11 different S*' for / > 1. 

In [4] , the three following states of polarization are studied. In the framework of the classical Stokes 
formalism, they are all considered to be identical states of polarization corresponding to totally depo- 
larized light. However, they all have different higher order statistics: i) A state of polarization with 
reduced Stokes vector S distributed uniformly on the Poincare sphere, ii) A state of polarization with 
reduced Stokes vector S distributed uniformly on the equator of the Poincare sphere (only linearly po- 
larized light), iii) A state of polarization with reduced Stokes vector S taking the value (0, 0, 1)* with 
probability 1/2 and the value (0, 0,-1)* with probability 1/2 (only left or right circularly polarized light). 

In all these three cases E(5) = — (0,0,0)* and P = P^ = 0. So that, in the classical Stokes 
formalism, they all correspond to the same state of polarization. Using the vectors 5', with I > 1, we 
can see how they are different: i) For this state, the vector is zero and P' = for all I > 1. ii) For 
this state 5'! = (0,0,0)* and P^ = 0, however = (0,0,-1/2,0,0)* and P^ = 1/2. iii) For this state 

51 = (0, 0, 0)* and Pi = 0, however 5^ = (0, 0, 1, 0, 0)* and P^ = 1. 

It appears from these three examples, in addition to the somewhat evident fact that higher order 
statistics are necessary when studying non Gaussian fields, that in order to call a state of polarization 
totally depolarized, it is not sufficient to have P = 0. Indeed, the state of polarization can verify P = 
but still refer to a specific type of polarization ellipse: only linear polarization (state ii)), or only circular 
polarization (state iii)). More examples are given in [4,5]. In [16] a new definition of the degree of 
polarization is proposed which takes into account higher order statistics. This definition is based on the 
following quantity: 



^ -ii^'ii-iiei 



(28) 



p' = \\e\\<vw+^)\a 

Where —l<n<lis such that maa;{„;<,„<;} [^^[ = [^^[. Now [10, 13]: 





(29) 
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Where p{s) is the probabihty density of the reduced stokes vector S. D is actually a quadratic measure 
of the difference between p(s) and a uniform distribution, whose probability density is equal to 1 . The 
degree of polarization is then defined as [16]: 



D 



e [0, 1] 



(30) 



P 



l + D 



This definition is closely related to the quantities P' , which we introducecd in the last subsection l5.ll In 
fact [10,13]: 



It is possible to synthesize the results of our model, with the measure of degree of polarization proposed 
in [16], by defining a totally depolarized state of polarization as one for which P' = for alH > 1. If 
we have = only for 1 < Z < L, then we can say that the state of polarization is depolarized to the 
order L. Classical models only consider depolarization to the order 1. Also, we should consider a state 
of polarization to be totally polarized only for P' = 1 for all I G N. This corresponds to a distribution 
concentrated at one point on the sphere. 

Let us make a final observation, without developing it: The measures of degree of polarization, and 
the criteria for distinguishing states of polarization with the same second order statistics, proposed in this 
article and in [4,5, 16] do not take into account the intensity distribution of the lightwave. The definition 
of degree of polarization proposed by Refregier [6,24], is based on the whole probability distribution of 
the field of the lightwave, including its intensity distribution. A comparative study of these two general 
approaches may help clarify the importance of including the intensity distribution in a measure of the 
degree of polarization. 

6 Evolution of the state of polarization during propagation 

In this section, we use our main equations, ([201) and of section 01 to study the evolution of the 
state of polarization of a lightwave, propagating in a random birefringent medium. This problem, as 
mentioned in the introduction, arises in optical fiber telecommunications [2,8, 17]. By using the model 
we have introduced in this paper, we will be able to achieve a detailed probabilistic description of the 
problem at hand. In particular, we will be able to give an exact analytical expression for the probability 
density on the Poincare sphere, representing the state of polarization, after any distance of propagation. 
We will also be able to follow exactly the evolution of different averaged quantities, such as the degree 
of polarization. 

The physical problem we are interested in is the one described by equation Q in section We have 
a lightwave, propagating in a random birefringent medium. We note z > the distance along which 
the wave has propagated. To each z > a random variable S{z) on the Poincare sphere is associated. 
It represents the state of polarization after a distance z in the medium. As explained in section |21 we 
suppose that there exists, for every z a rotation R{z) such that: 



The evolution of the state of polarization during propagation can then be described by a stochastic 
process S{z) on the sphere, or a stochastic process R{z) on SO{3). Our description of this evolution is 
based on the mathematical concept of a Levy Process on S0(3) [14], which we introduce in the following 
subsection 16.11 

6.1 Levy process model 

We model R{z) as a Levy process on 5*0(3). This model reflects a set of simple physical properties of 
the propagation medium. It is based on the following hypotheses: 

• Independent increments: For zi < Z2 we have that R{zi) and R{z2)R^^{zi) are independent. 
Physically, this means that non overlapping parts of the medium are not coupled. 




(31) 



i>i 



S{z) = R{z)S{0) 



(32) 
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• Stationary increments: For zi < Z2 we have that R{z2)R^^{zi) = R{z2 — zi). Physically, this 
means that the medium is homogenous and only locally random. This hypothesis, more generally 
means that: 

R{z) = R{z - Zn)R{Zn - Zn-l)...R{z2 - Zi)R[zx) for Z > Z„ > Z„_i > ... > Z2> Zy (33) 

• Stochastic continuity: The stochastic process Riz) is stochastically continuous. This means 
that the probability for R{z\) and R{z2) to be different tends to zero as Z2 — z\ goes to zero. 
Physically, this means that a very short length of the medium can not induce a big change in the 
state of polarization^ 

• We add the simplifying hypothesis that i?(0) = /, where / is the 3x3 identity matrix. 

A stochastic process on 50(3) -or any other Lie group- verifying these properties is called a left Levy 
process, or just a Levy process [14]. 

This model can be very effectively reduced, using the main equations ([20|) and ([23ll . Note p_R.(z)(r), 
the probabihty density of R{z^. According to equation 



R{z) = R{z - zi)R{zi) for z > zi (34) 
Using the results of section U] -equation (PTjl - we can write: 

PB.(.z)ir) = iPB^{z-z^) *PR(z,)){'r) (35) 

Using equation ()22p . we transform this last equation: 

R\z) = R\z - zi)R\zi) (36) 

Where i?'(z) are the coefficient matrices in the development of Pr[z)- 

The stochastic continuity of the process R{z) implies the continuity in z of the matrices R\z). The 
only continuous solution of ([55]) verifying i?(0) = / is [14]: 

i?'(z) = e*'Mor / > 1 ^ ' 

Where (for every Z > 1) is a constant matrix (not function of z). These constant matrices are called 
generators of the process R{z). We have: 

= ^R'i^)\z=o (38) 



It follows from the decomposition formula p4p that: 

PRiz)ir) = l + J2i2l + l)Trie''^iD\r)y) (39) 
i>i 

Which gives the probability density of R(z) fr any z. 

From this last result, the probability density on S'^ of S{z) can be derived in a direct way. We have 
already noted that S{z) = R{z)S{0). H Ps{z){s) is the probability density of S{z) and f'(z) are the 
cefficient vectors in its development (fT^ then by equations (fT^ and ((^ : 



C"(z) = 1 

^'(z) = e*"^^'(0) for / > 1 



(40) 



Using the decomposition formula (|T3|) it follows that: 

Ps(z){s) = 1 + + l)(e*'^^'(0))*F'(s) (41) 



i>i 



^This does not impose that the trajectories on the Poincare sphere of s(z) are continuous. They can have jump 
discontinuities. 
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Which gives the probabihty density of S{z) for any z. 

By modelhng the evolution of the state of polarization, during propagation in a random medium, as 
a Levy process, the descriptions (f5^ and ([iT|) of this evolution have been achieved. Practically, these 
descriptions give the evolution of the state of polarization in function of the generator matrices t', > 1). 
It is clear from equation ([38|) that these matrices characterize the propagation medium locally. That is, 
they describe propagation over small " ifinitesimal" distances. These matrices are not known a priori. 
There are two ways of giving them: The first way is to use a local model for the propagation, in the form 
of a stochastic differential equation. Such local models can be found in [2,17]. The second way is to note 
that these matrices are parameters of the probability density of Ps(z) ■ It is possible to estimate them 
given realizations of S{z) -see discussion in the introduction. We will return to these two approaches in 
subsection 16.31 

6.2 Depolarization 

The most important effect of a random medium on the state of polarization of a lightwave is to depolarize 
it. After a long distance of propagation, we can expect the state of polarization of the lightwave to become 
totally depolarized. Here, we use the mathematical description presented in the last subsection 16.11 to 
give the evolution, during propagation in a random birefringent medium, of the degree of polarization. 
We see that this evolution tends to a totally depolarized state independently of initial conditions. 

It is possible to show, under very general conditions [14], that the (real parts of the) eigenvalues of 
the generator matrices in equation ((57)) are all negative. Since = e*'^, this means that i?' ^ 
exponentially for large z. By taking the limit of equation ([^T|) for large z, we find that Ps{z){s) = 1 
for large z. In other words the probability distribution of S{z) tends to a uniform distribution on the 
Poincare sphere, which is characteristic of a totally depolarized state. Note, from equation (|41|1. that 
as the matrices tend to zero Ps{z){s) will tend to a uniform distribution independently of the initial 
distribution Ps{o)is)- 

We have established that the state of polarization of a lightwave propagating in a random medium 
tends, with the distance of propagation, to a totally depolarized state represented by a uniform probability 
distribution on the Poincare sphere. Let us now examin the dependence on z of the degree of polarization 
of order P', for Z > 1. Remember that was defined in section [?TT1 equation (pS)) . as the norm of the 
complex vector Using equation (j40p . we can write: 

p\z) = \\eiz)\\ = \\e''^em m 

It is clear in this equation that P^{z) tends to zero independently of initial conditions. This is, in 
particular, true for the usual degree of polarization P = P^. 

We have just given the law of evolution, i.e. the dependence on z, of the degree of polarization of 
order I, P'', for all / > 1. These quantities are an example of what might be called averaged quantities 
associated to the process S{z). That is, combinations of the averages of functions of S{z). Averaged 
quantified are, of course, deterministic. Let f{s) be any real square integrable function on the sphere 
S^. An example of an averaged quantity is E(/(S'(z))). The fact that / is real can be used to slightly 
transform its development (jl3p . in the following way: 

/(.) = M = ^(2/ + myy^) m 

Using the fact that, equation £,''{z) = E (Y^{s)j , it is possible to write: 

E(/(5(z))) = 5](2/ + l)(/')te'(z) (44) 

Or, using the expression PO)) for ^'(z): 

E (/(5(z))) = 5^(2/ + l)ih^e''^e{0) (45) 

The last equation (|45l shows that the evolution of any averaged quantity can be followed exactly if the 
generator matrices arc known. Averaged quantities include the entropy of the state of poalrization. 
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the average parameters of the ehipse of polarization, or any other attribute of the state of polarization 
that we may wish to study. Equation (|45p can for instance be used to establish that the informational 
entropy of the state of polarization is stricly increasing during propagation in a random birefringent 
medium. Thus, depolarization can be associated with an increasing informational entropy. 



6.3 Estimation of physical parameters: An example 

In subsection 16. 11 the evolution of S{z) during propagation was modelled as a Levy process. This model 
lead to an analytical formula (j4ip . for the probability density of S{z), containing the generator matrices 
t' {I > 1) - see equations (|57)) and (IHS)) - as free parameters. As mentioned before, these matrices are 
related to the local properties of the propagation medium. In order to give them concrete expressions or 
values, two approaches can be used: The first is to use a local physical model for the evolution of S{z). 
This model can take the form of a stochastic differential equation [2, 17]. The second way is to estimate 
them from realizations of S{z), since indeed, these matrices appeare as parameters of the probability 
density of S{z). 

A particularly simple case arises when the generator matrices correspond to the stochastic equation 
proposed in [2]. This is a stochastic differential equation describing the evolution of the vector on the 
Poincare sphere S{z) during propagation in an optical fibre affected by PMD. The main idea of this 
model is that S{z) rotates on the Poincare sphere with an angular velocity which is essentially a white 
noise vector: 

^S'(z) = fiW X Siz) (46) 

Where /i is a constant, is a three dimensional white noise vector, x is the vector (cross) prodcut and 
the equation is to be understood as a Stratonovich stochastic differential equation [2] . 

This stochastic differential equation is well known in mathematics [7,13,14,21]. It describes Brownian 
motion on the sphere . The generator matrices for this process are given by [7, 13, 14,21]: 

t' = ^/(/ + l)/, (47) 

Where /; is the {21 + 1) x (21 + 1) identity matrix. By replacing this result in (|4T|) . it follows that: 

Ps(.)(s) - + l)e"^'('+^'^(^'(0))*r'(s) (48) 

This situation is particularly simple. By starting from the local model suggested in [2], we arrive at 
an expression for the probability density of S{z) which depends on only one free parameter, namely /Lt, 
instead of having the (infinite) family of matrices t' as free parameters. 

note also that, in this case, the degree of polarization of order / takes on a simpler form: 

p'(z) = e^'('+i)^||C'(0)|| (49) 

Now let us see how, in this simple case which only has one parameter in the probability density of >5'(2:), 
we can use realizations of S{z) to estimate /i. We need to consider the medium as a whole (a closed 
system). We note S{0) = Sin and consider a length Z of the medium so that we can note S{Z) — Sout- 
The medium is represented by the rotation R = R{Z), which is the essential part of its Mueller matrix. 
In an experimental framework Si„ should be known to us. Here, it is assumed to be a pure state of 
polarization such that Sin takes the value (0,0, 1)* with probability one (left circular polarization). In 
this casell, the probability density of Sout is given by formula fiSl) : 

= + l)e^'('+i)^yo'(s) (50) 

Where we have replaced the values of C'(0) corresonding to Sin- And the degree of polarization is given 
by- using formula (|49|l : 

P = Fi = e-^'^ (51) 



'This formula was first given by F. Porrin in 1928. 
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Formula (jSOp gives the probability density function of Sout ■ This probability density function contains /i^ 
as a parameter. A standard way of estimating /i^ is, for instance, maximum of liklihood estimation [9]. 
Using formula (|50p in order to find an analytical expression of the maximum of liklihood estimator of 
/i^ is not a straightforward task. Still, an exact maximum of liklihood estimator of /i^ can be found 
numerically in a standard way. When ^^Z is small, the maximum of liklihood estimator of /i^ can be 
approximated with the following estimator [15,20]: 



i=N 

2ZN ^ 

i=l 



(52) 



The estimator is evaluated from realizations Si -with i = 1, 2, . . . , A^- of the random variable Sout- Here 
9i is the polar angle of 5*^. That is, its angular distance, on the Poincare sphere, from the inial value 
Sin — (0,0,1)*. The estimator is thus based on the emperical mean of the squared angular distance 
between the initial state Sin and each realization of Sout- An estimator similar to this one is used in [18] 
for the standard deviation of speckle noise. Another approach to the estimation of fi^ is to estimate the 
degree of polarization and use equation ([ST]) to retreive /i^ . According to our definition of the degree of 
polarization, formula (I25p . P can be estimated as: 



i=N 



(53) 



The degree of polarization can also be estimated from intensity measurements [23] 



7 Conclusions and outlook 

This article was aimed at presenting a new model for the propagation of polarized light in random 
birefringent media. The physical situations to which this model would be applicable are bounded by 
the hypotheses introduced in section [5] This model was intended to be adapted to a detailed statistical 
treatment of the physical problems it describes. This is done by including higher order statistics of the 
reduced Stokes vector and by describing the random medium via a certain type of statistical input/output 
relation. It was argued that this type of relation, as introduced in sectionHl can accomodate a variety of 
physical models and make it easier to estimate the physical parameters appearing in these models. The 
model that was presented is based on a decomposition of the higher order statistics of the reduced Stokes 
vector along the irreducible representations of the rotation group 5*0(3), which is the group giving the 
action of a birefringent medium on the Poincare sphere. In section 15.11 this decomposition was used 
to generalize the notions of reduced Stokes vector and degree of polarization to higher order statistics. 
This generalization was discussed in relation to recent works studying the role of higher order statistics 
in polarization optics. In addition to this more theoretical result, the decomposition along irreducible 
representation was used, section^ to give the laws of transformation, i.e. the input/output relations, for 
the higher order statistics of the reduced Stokes vector of a lightwave propagated through a birefringent 
random medium. These laws of transformation arise mathematically from the spherical convolution 
theorem. 

In section[6l the evolution of the state of polarization of a lightwave propagating in a random birefrin- 
gent medium was studied. It was modelled using the concept of Levy Processes on the rotation group, 
see subsection 16. II The framework of Levy processes on the rotation group was used to give an analytical 
expression of the probability density on the Poincare sphere, representing the state of polarization after 
propagation over any distance in the medium. In subsection 16.2) the Levy process model was used to 
study the depolarization of a lightwave by propagation in a random birefringent medium. An analytic 
law for the evolution of the degree of polarization under the effect of propagation in the medium was 
given. It was shown that depolarization takes place independently of initial conditions. In subsection 
16.31 an example of propagation in optical fibres was discussed. In particular, the probability density on 
the Poincare sphere, representing the state of polarization of a lightwave propagating in an optical fibre 
affected by PMD, was given in an analytical form. The estimation of the physical parameters appearing 
in this prbability law, from observations of the reduced Stokes vector was discussed. 

An important question to adress is how the model presented in this paper, which is specific to 
birefringent random media, can be generalized to any random medium. Such a general random medium 
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acts on the Poincare sphere by nonhnear transformations. This makes it difficult to model using the 

Poincarc sphere formalism. However, if wc use the complete Stokes formalism (i.e. With all the four 
components of the Stokes vector), then the medium acts on the Stokes vector essentially by Lorentz 
transformations [1,22]. These are, of course, linear transformations. 

In this article, we have deduced the laws of transformation of the higher order statistics of the 
reduced Stokes vector by using the irreducible representations of the rotation group. In the case of a 
general random medium, it is possible to do the same for the higher order statistics of the complete 
Stokes vector. The finite dimensional irreducible representations of the Lorentz group (which contains 
the rotation group as a subgroup) should then be used. The finite dimensional representations of the 
Lorentz group are known as spinor representations [19]. All the main features of the model presented in 
this paper can be, in this way, generalized to any random medium. However, somme additional technical 
difficulties would arise, since the Lorentz group, unlike the rotation group, is not compact. 

The general idea of the model we have presented is to group as much a priori knowledge as possible, 
about polarized light in random birefringent media, in one consistent statistical framework taking into 
account the higher order statistics of the wave field. In other words, to construct a general signal model 
for the state of polarization of a lightwave in a random birefringent medium. This signal model would 
allow the extraction of significant physical information, in a variety of practical problems, using different 
signal processing techniques such as detection, estimation, filtering, etc. In this paper, we have presented 
our signal model and given a toy example of how it can be used. We hope to demonstrate the usefuleness 
of this model by applying it to concrete problems in future works. 
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